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Kj^.^! Abstract: Bolometers have proven to be good detectors for the search of neutrinoless double beta 

J^ \ decay. By operating at cryogenic temperatures, they feature excellent energy resolution and low 

■ ' background. The detection of the possible light emitted when particles interact in the bolometer 

is a promising method to lower the background of the experiments. The different amount of light 

^ ' emitted in j3/7 and a interactions, whether due to scintillation or Cerenkov emission, allows to 

^ ■ discriminate the two interaction types. Because of the cryogenic environment, light detectors are 

\^ . often bolometers. In this work we present a software algorithm to lower the energy threshold of 

bolometric light detectors coupled to luminescent bolometers. The application to data from Ge 

light detectors coupled to ZnMo04 and Te02 bolometers shows that the energy threshold can be 

lowered substantially, increasing the discrimination power when the amount of emitted light is 

small. 
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1. Introduction 

Bolometers are detectors in which the energy from particle interactions is converted to heat and 
measured via their rise in temperature. They provide excellent energy resolution, though their 
response is slow compared to electronic or photonic detectors. They are used in particle physics 
experiments searching for rare processes, such as neutrinoless double beta decay (OvDBD) and 
dark matter (DM) interactions. 

The CUORE experiment will search for OvDBD of i^°Te % g] using an array of 988 Te02 
bolometers of 750 g each. Operated at a temperature of about 10 mK, these detectors maintain an 
energy resolution of a few keV over their energy range, extending from a few keV up to sev- 
eral MeV. The measured resolution at the Q- value of the decay (Q = 2528 keV [^]) is about 
5 keVFWHM; this, together with the low background and the high mass of the experiment, de- 
termines the sensitivity to the OvDBD. 

To further increase the sensitivity, an intense R&D is being pursued to lower the background 
in the OvDBD region. The main source of background is due to a particles, coming from radioac- 
tive contaminations of the materials facing the bolometers. The natural way to discriminate this 
background is to use a scintillating bolometer |Q]. In such a device the simultaneous and indepen- 
dent read out of the heat and the scintillation light permits to discriminate events due to jS/y and 
a interactions thanks to their different scintillation yield. Unfortunately Te02 is not a scintilla- 
tor, and different crystals are being studied. Some of them are CdW04 (OvDBD candidate '^^Cd, 
Q = 2902 keV) [|], ZnSe (^^''Se, Q = 2995 keV) [|], and ZnMo04 ("'°Mo, Q = 3034 keV) [0]. 

To detect light the scintillating crystal is faced to a high sensitivity dark bolometer ^. It usu- 
ally consists in a germanium or silicon slab that absorbs the scintillation photons giving rise to a 
measurable increase of its temperature. The energy threshold of these devices plays an important 



role in the identification of a and j8/7 interactions. ZnMo04, for example, is a promising candi- 
date [^ even though the scintillation yield is one of the smallest among the tested crystals. In a 
30 g ZnMo04 bolometer, in fact, only about 1 keV of light is emitted per MeV of particle energy 
by P/y's, and about a tenth by a's. Currently the energy threshold of light detectors is around 
300 eV, so that a interactions cannot be completely studied at 3 MeV, where the OvDBD signal is 
expected. Moreover increasing the mass of the bolometer to the kg scale, as needed by a CUORE 
like experiment, could substantially reduce the scintillation yield, reducing the separation between 
P/Y's and a's. 

In this paper we present a software algorithm that lowers the threshold of the light detectors. 
The application to data from a ZnMo04 bolometer shows substantial improvements. The algo- 
rithm was also applied to a TeOi bolometer coupled to a bolometric light detector. Thanks to this 
work, Ught, probably due to Cerenkov emission, has been detected for the first time, opening the 



possibility of rejecting the a background in a Te02 experiment [10|. 



2. Experimental setup 

The data presented in this paper come from a ZnMo04 crystal and a Te02 one, operated as 
bolometers in the Gran Sasso National Laboratories in Italy (LNGS), in the CUORE R&D facility 



TheZnMo04 crystal weighed 29.9 g and was a parallelepiped of dimensions 29 x 18 x 13 mm^. 
The light detector (LD) consisted of a 36 mm diameter 1 mm thick pure Ge crystal, covered with 
a 600 A layer of Si02 to ensure good light absorption. A reflecting foil (3M VM2002) was placed 
around the ZnMo04 crystal to enhance the light collection. 

The Te02 crystal weighed 1 16.7 g and was a parallelepiped of dimensions 30 x 24 x 28 mm^. 
The crystal was doped with natural samarium (see details in Ref. [|T^), which contains 15% of 
^"^^Sm, a long living isotope (T1/2 = 1.07 x 10^' y [p3|]) which undergoes a decay with a Q-value 



of 2310±1 keV. This decay allows a direct analysis of the behavior of a's in an energy region close 
to the OvDBD. As for the ZnMo04, the crystal was surrounded by a reflecting foil. The light was 
detected with a pure Ge bolometer larger than the one faced to the ZnMo04, its dimensions being 
66 mm diameter and 1 mm thickness. 

The temperature sensor of the ZnMo04 crystal was a Neutron Transmutation Doped (NTD) 



germanium thermistors [16, 17] of 3x1.5x0.4 mm^, thermally coupled to the crystal surface by 
means of 6 epoxy glue spots of about 0.6 mm diameter and 50 /im height. For the Te02 the 
NTD-Ge thermistor had dimensions 3x3x1 mm^, and was glued with 9 glue spots. Each LD 
was equipped with two NTD-Ge thermistors of 3 x 1.5 xO.4 mm^. The ZnMo04 and Te02 crystals 
and the LDs were held in a copper structure by Teflon (PTFE) supports, thermally coupled to the 
mixing chamber of a dilution refrigerator which kept the system at a temperature of about 13 mK. 
The read-out of the thermistors was performed via a cold pre-amplifier stage, located inside 
the cryostat, and a second amplification stage, located on the top of the cryostat at room temper- 
ature. After the second stage, the signals were filtered by means of an anti-aliasing 6-pole active 
Bessel filter (120 db/decade), and then fed into a NI PXI-6284 analog-to-digital converter (ADC) 
operating at a sampling frequency of 2 kHz. 
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The trigger was software generated on each bolometer. When it fired, waveforms 0.6 s long 
on the ZnMo04 and 2 s long on Te02 were saved on disk. The length of the window acquired 
from the LDs was equal to that of the crystal they faced. If the trigger fired on the ZnMo04 or on 
the Te02, waveforms from the corresponding LD were anyhow saved, irrespective of their trigger. 
To maximize the signal to noise ratio, waveforms were processed offline with the optimum filter 



algorithm [18, 19] 



The main parameters of the bolometers are reported in Tab. [l[ The rise and decay times of 
the pulses are computed as the time difference between the 10% and the 90% of the leading edge, 
and the time difference between the 90% and 30% of the trailing edge, respectively. The intrinsic 
energy resolution of the detector is evaluated from the fluctuation of the detector noise, after the 
application of the optimum filter. 

Table 1. Parameters of the bolometers. Amplitude of the signal before amplification (As), intrinsic energy 
resolution after the application of the optimum filter (a), rise (Tr) and decay (t^) times of the pulses (see 
text). 





As 
[AtV/MeV] 


a 

[keV RMS] 


T, 
[ms] 


[ms] 


ZnMo04 

Ge (ZnMo04) 


321 
1.7x10^ 


0.5 
0.07 


8 
4 


33 
11 


Te02 

Ge (Te02) 


43 
2.3 xlO^ 


1.31 
0.10 


15 
3 


116 
10 



To perform the energy calibration and to provide statistics to the P/y sample the detectors 
were exposed to a ^^^Th source. The source generates /'s with energy up to 2615 keV, the energy 
of a y-ray from ^°^T1. The Te02, together with the 2.3 MeV a line from ^"*^Sm, also featured 
a contamination in ^^"Po, which undergoes a decay with a Q-value of 5.4 MeV. To evaluate the 
discrimination between jS/y's and a's over a wide energy range, a smeared a-source was faced to 
the ZnMo04 crystal. The source consisted in 1 /il of a 0.1% ^^^U liquid solution, covered with a 
6 /xm mylar foil to absorb part of the a's energy and produce a continuum spectrum in the range 
1-4 MeV. As in the Te02 case, an a population at 5.4 MeV from ^^"Po contamination is visible in 
the ZnMo04 energy spectrum. The LDs were permanently exposed to a ^^Fe source, placed on the 
LD surface opposite to the crystal. The source produces two X-rays at 5.9 and 6.5 keV, which are 
used for the absolute calibration of the bolometer. 

3. Data analysis 

As said in the previous section, waveforms acquired by the ADC are processed with the optimum 
filter. This filter maximizes the signal to noise ratio and improves the energy resolution and the 
energy threshold of this kind of bolometric detectors (see for example Refs. [ pO| , pTj]). The transfer 
function is built using the signal shape s{t) and the noise power spectrum N{o)): 

H{w) = h'-^^e-J-'-, (3.1) 
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Figure 1. Waveforms detected on the ZnMo04 (top) and light (bottom) detectors for a 2615 keV 7-ray (left) 
and a 200 keV energy release (right). Original (dark line) and optimum filtered (red line) waveforms. The 
dashed line marks the maximum position of the filtered signal on the ZnMo04. 



where his a normalization constant, s{(o) is the Fourier transform of s{t) and tM is a parameter to 
adjust the delay of the filter. The signal shape is estimated by averaging a large number of triggered 
pulses. On the ZnMo04 and TeOi bolometers the pulses are selected in the high range of the 
energy spectrum (Energy > 1 MeV), on the LDs the pulses are selected from the peaks of the ^^Fe 
source. The noise power spectrum is estimated by averaging the power spectra of a large number 
of waveforms not containing pulses. 

Once the waveforms are filtered, the amplitude of a signal is estimated as the maximum of the 
filtered pulse. As said before, when the trigger fires on a heat bolometer, the waveform from the 
corresponding LD is acquired. The amplitude of the light signal is estimated, also in this case, as 
the maximum of the filtered waveform. Signal amplitudes are then corrected for temperature and 



gain instabilities of the setup [ ]22[ ] and calibrated in energy. 

Figure |l] shows two pulses triggered on the ZnMo04 bolometer and the corresponding wave- 
forms acquired from the light detector. The pulse on the left is generated by a 2615 keV 7-ray 
interaction, which also induces a visible pulse on the light detector. The pulse on the right is gen- 
erated by an interaction of a 200 keV particle (yor j8) which, given the low amount of scintillation 
Ught produced at this energy, does not induce a clear pulse on the light detector. While in the first 
case the maximum of the filtered waveform results in a correct estimation of the light signal ampli- 
tude, in the second case the maximum could be a fluctuation of the noise, rather than the amplitude 
of the signal. Therefore, when the amount of emitted light is at the level of the noise or below, the 
estimation of the signal amplitude with a maximum search algorithm may produce wrong results. 

This effect is visible in the distribution of the light detected versus particle energy (Fig. ^). At 
high energies the light emitted from 7 interactions is proportional to the particle energy. At low 
energies, below 300 keV, the detected light reaches a constant value at about 250 eV, which is the 
noise pedestal. This effect is also visible at higher energies, in the interaction of a particles which 
produces less scintillation light than jS/yones. In this case the pedestal is reached when the particle 
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Figure 2. Light detected from particle interactions in the ZnMo04 crystal using the maximum search 
method. 



energy is below 3 MeV. As it will be shown in Sec. ^^ the Te02 case is even more dramatic. The 
amount of emitted light is so small that the noise pedestal is measured for all particle energies, 
irrespective of their jS/y or a nature. 



4. Lowering the energy threshold 

The bolometers used in this work are detectors that respond in milliseconds, a time orders of mag- 
nitude larger than the process of light emission. Therefore the energy releases into the heat (HD) 
and light (LD) detectors can be considered synchronous. The time delay between the two signals 
depends only on the differences between the thermal responses and the readout circuits of the two 
bolometers, and hence is fixed. 

The basic idea of this work is to estimate the amplitude of the signal contained in the LD 
waveform as the value of the filtered waveform at a fixed time delay with respect to the position of 
the maximum in the HD waveform. 

The time delay can be estimated from a set of events in which the energy released on both 
bolometers is much higher than the noise. In this case the maximum search algorithm does not 
fail on the LD, and the time difference between the maximum positions of the signals on the two 
bolometers is a good estimator of the delay. Taking the ZnMo04 as example, the time delay is 
estimated from events generated by high energy jS/y's, like the one in Fig. |l| left. Then, instead of 
looking for a maximum in the LD filtered waveform, the amplitude value found at the estimated 
time delay with respect to the maximum position in the HD is used. 

4.1 Application to the ZnMo04 bolometer 

The time delay is estimated selecting the events in Fig. ^ with energy in the range 1000-2650 keV 
and light energy in the range 1.0-2.5 keV. The distribution of the time difference between the 
maximum positions of the signals in the ZnMo04 and in the LD is shown in Fig. ^ The time delay 
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Figure 3. Distribution of the delay between the light signal and the ZnMo04 heat signal using the events in 
the yband in Fig. |[ The events selected lie in the energy range 1000-2650 keV and in the light energy range 
1.0-2.5 keV. The plot shows all the events found in the allowed range of time differences. The inset shows 
a zoom of the peak. 



is not known a priori, so a wide range of time intervals is allowed in the software algorithm. The 
average of the distribution is at -9 ms, a value which is expected to be negative, since LDs responds 
faster than HDs (see the rise times of the bolometers in Tab. |l|). 

In the wide range of times allowed, spurious events lying off the distribution core may occur. 
We chose to estimate the time delay as the mode of the distribution, which is an estimator resilient 
to outliers. The average or a Gaussian fit, in fact, are not resilient to outliers, and may produce 
incorrect results when implemented in automatic software routines. In the ZnMo04 case the mode 
coincides with the average, because there is only one outlier. As it will be shown in the next section, 
the outliers in the Te02 distribution are many and spoil the distribution. 

We applied our method to the same waveforms used to produce the scatter plot in Fig. ^, 
and we show the comparison in Fig. ^. While at high jS/y energies the two methods produces 
compatible results, when the energy goes below 300 keV the light energy estimated with the new 
method goes below the noise pedestal and is statistically zero for zero energy released in the HD. 
The same behaviour is visible on a particles. 



Bolometers are nonlinear detectors [|23[], implying that the shape of the signal depends on the 
energy released. The time delay depends on both the energies released in the HD and in the LD, 
and in principle the same value cannot be applied to the entire range of energies. We checked this 
and we observed that the variation of the time delay is small and produces negligible changes in 
the estimation of the LD amplitude. 

The reader should have noticed that the light energy estimated with the new method can take 
negative values. At a first glance this may seem odd, but it is not. Bolometers are thermal detectors. 
Fluctuations of the detector baseline corresponds to temperature variations, which may by positive 
or negative. When the signal is smaller than the noise or null, the estimated energy can be positive 
or negative as well. 

To quantify the improvements reached with the application of our method to the ZnMo04 
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Figure 4. Light detected from the ZnMo04 crystal using the maximum search method (black dots) and the 
new method (red dots). 

bolometer, we study the threshold of the LD using the events produced by very low energy jS/y 
particles (events with Energy < 20 keV in Fig. ^. In these events the light emitted is negligible, 
and the energy measured by the LD is dominated by the noise, which sets the ultimate energy 
threshold. Erom the distributions shown in Eig. ^ we observe that the light energy measured with 
the maximum search algorithm has an asymmetric shape, while the energy measured with the new 
method has a Gaussian shape. The a of the noise Gaussian, estimated by means of a fit, is found 
to be compatible with the intrinsic resolution of the light detector. The 50% of the noise events is 
below 249 eV using the old method and below 19 eV using the new method. The 90% is below 
382 eV using the old method and below 125 eV using the new method. 

4.2 Application to the TeOi bolometer 

The time delay cannot be estimated for the TeOi as in the case of ZnMo04, since, as shown in 
Eig. |6| (top), the maximum search algorithm clearly fails to find the signal in the LD. The noise 
pedestal is measured for all particle energies. 

Lacking the signal sample, we estimated the time delay using events generated by particles 
interacting in both the LD and the HD. Also in this case, even if the nature of the energy release 
in the LD is different from light, the energy releases on the two bolometers can be considered 
synchronous. To avoid accidental coincidences, and eliminate noise fluctuations and spikes, we 
selected events with energy greater than 200 keV on both detectors. The events selected, due to 
multiple interacting y's, and the distribution of their time delay are shown in Eig. ^. The outliers in 
the time delay distribution are due to residual accidental coincidences between particles interacting 
in the detectors. The mode is at -27 ms which is assumed as the value of the time delay. We checked 
on the ZnMo04 that the time delay estimated with events generated by particles interacting on both 
detectors is compatible with the delay of the light signal. 

Erom the distribution of the light estimated with the new method versus particle energy (Eig. ^ 
bottom) we see that a clear correlation appears in the jS/y band. This allows to discriminate these 
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Figure 5. Distribution of the light detected from the ZnMo04 crystal applying the selection Energy < 20 keV 
to the events in Fig. H, using the maximum search method (top) and the new method (bottom). The blue solid 
and red dashed lines represents the 50% and 90% of the distributions integrated from the left, respectively. 
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Figure 6. Light detected from the TeOa crystal using the maximum search method (top) and the new method 
(bottom). The events above the 2615 keV line are a particles, the others are j3/7 particles, except for the 
line at 2310 keV (green dots) which is an a decay from ''^^Sm in the Te02 crystal. 



particles from a's which do not produce a detectable energy. We estimate the threshold of the LD 
using low energy events (events in Fig. ^ with energy less than 200 keV). The 50% of the noise 
events is below 328 eV using the old method and below 4 eV using the new method. The 90% is 
below 437 eV using the old method and below 144 eV using the new method (Fig. P). 
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Figure 7. Distribution of the time delay between the responses of the Te02 crystal and the light detector 
(bottom left) on events generated by particles interacting in both detectors. The events are selected requiring 
that the energy released is greater than 200 keV on both detectors (top). The mode of the delay distribution 
(right) is -27 ms. 
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Figure 8. Distribution of the light detected from the TeOi crystal applying the selection Energy < 200 keV 
to the events in Fig. ^ using the maximum search method (top) and the new method (bottom). The blue solid 
and red dashed lines represents the 50% and 90% of the distributions integrated from the left, respectively. 



5. Conclusions 

The algorithm we developed improves considerably the results that can be obtained from lumi- 
nescent bolometers. The search of the light signal through the maximum of the optimum filtered 
waveform produces wrong results when the signal is at the level of the noise or below. In place of 
the signal amplitude, the noise pedestal, which is higher, is measured. The new method uses the 
heat bolometer as reference. The amplitude of the light signal is estimated from the value of the 
filtered waveform at a fixed time delay with respect to the signal in the heat bolometer. In this way 
the noise pedestal is eliminated. The amplitude measured is the real signal amplitude, smeared by 
the noise distribution, which is a Gaussian centered at zero amplitude. With the new method the 
energy threshold corresponding to a 10% rate of noise resembling signal is reduced by a factor 
about 3. 

Finally, an important feature of our algorithm is that, even if the signal is below the noise, 
it produces an unbiased estimation of the signal amplitude. By increasing the statistics any small 
amount of light can be estimated, which is useful to study the features of the light emitted by 
particles of different nature and energy. 
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